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A new, improved version of a cosmic crystallography method for constraining cosmic topology is introduced. Uike the circles-in- 
the-sky method using CMB data, we work in a thin, shell-like region containing plenty of objects. Two pairs of objects (quadruplet) 
linked by a holonomy show a specific distribution pattern, and three filters of separation, vectorial condition, and lifetime of objects 
extract these quadruplets. Each object P, is assigned an integer .s,, which is the number of candidate quadruplets including P, as their 
members. Then an additional device of i,-histogram is used to extract topological ghosts, which tend to have high values of Sj. In 
this paper we consider flat spaces with Euclidean geometry, and the filters are designed to constrain their holonomies. As the second 
filter, we prepared five types that are specialized for constraining specific holonomies: one for translation, one for half-turn corkscrew 
motion and glide reflection, and three for ;i-th turn corkscrew motion for n = 4, 3, and 6. Every multiconnected space has holonomies 
that are detected by at least one of these five filters. Our method is applied to the catalogs of toy quasars in flat A-CDM universes 
whose typical sizes correspond to z ~ 5. With these simulations our method is found to work quite well. These are the situations 
in which type-II pair crystallography methods are insensitive because of the tiny number of ghosts. Moreover, in the flat cases, our 
method should be more sensitive than the type-I pair (or, in general, n-tuplet) methods because of its multifilter construction and its 
independence from n. 
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1. Introduction 

The shape of our space is one of the greatest and the oldest ques- 
tions in human history. Ancient people tried to answer the ques- 
tion mythologically or philosophically, but in the present era, we 
are ready to answer it scientifically. When we consider the uni- 
verse as a 4-manifold consisting of 3-space and 1-time, answer- 
ing the question is translated into determining the local geometry 
and the global topology of our Universe. 

Local geometry is described by Einstein's General Relativity 
(GR). The metric of space-time plays a fundamental role, and the 
assumption of the cosmological principle, which states that our 
Universe is (locally) homogeneous and isotropic, leads us to a 
famous Friedman-Lemaitre-Robertson- Walker (FLRW) metric, 

ds^ = -c^dt^ + R{tf\ + r^ide^ + sin' ed4>^)\ (1) 

L 1 - kr^ J 

where k is the normalized curvature of space, such as A; - +\ 
(spherical geometry), (Euclidean geometry), or -1 (hyperbolic 
geometry). Together with Einstein's equations, curvature is re- 
lated to the average energy density of the universe and then to 
other physical quantities. Recent observations favor a A-CDM 
universe with curvature k ^ Q (e.g. Qtot - 1 .0050;^Q{*gj from 
WMAP+'&AD+SN data, by Hinshaw et al. 2009), suggesting that 
our Universe is one of the flat spaces with Euclidean geometry. 

Global topology, on the other hand, has no reliable physical 
theories to describe it, so the easiest approach is to constrain 
it mathematically through direct observations. This situation is 
analogous to the case of Carl Friedrich Gauss who had to mea- 
sure the angles of the large triangle formed by three peaks of 
mountains to know which geometry describes our space, since 
he did not know GR and modern cosmology. 



Moreover, it is true that the global topology of the universe 
has been a relatively less popular concept than the local ge- 
ometry. A pioneering work on cosmic topology was done by 
Ellis (1971). Other early works include Sokolov and Shvartsman 
(1974), Fang and Sato (1985), Gott (1980), Fangundes (1983), 
and so on. This theoretical or observational research has not at- 
tracted much notice due to the lack of observational data. In the 
past two decades, however, we have seen tremendous progress 
in this field, along with progress in observational techniques. 
Specifically, a possibility of the Poincare dodecahedral space 
topology suggested by Luminet et al. (2003) was a breakthrough 
in cosmic topology. This field recieves more and more interest 
these days, not only from theorists but also from observational 
astronomers. The overall topology of the Universe is now be- 
coming one of the major concerns in astronomy and cosmology. 

In modem cosmology with an FLRW metric, already de- 
scribed above, we have three geometries (scalar curvatures) ac- 
cording to their sign. For each geometry, there is only one type 
of space with a simply connected topology, namely, 3-sphere S^, 
3-Euclidean space E^, and 3-hyperbolic space H^, for k - +1,0, 
and -1, respectively. It is possible to construct a space that is lo- 
cally indistinguishable from the simply connected one, i.e., hav- 
ing the same curvature but with a different topology. It is not 
simply connected, but multiconnected. 

We give a brief review of multiconnected spaces below. 
Detailed treatments are found in various reviews (e.g. Lachieze- 
Rey and Luminet 1995). A multiconnected space M is a quotient 
space of the simply connected space U with the same geome- 
try, by a holonomy group F, where a holonomy is an isometry 
on U without any fixed points (except for the identity). Hence, 
M - U/r can be imagined as U tiled by polyhedra identified 
by holonomies y e F, possessing repeated structures. This poly- 
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hedron, called a fundamental cell, is a 2/r-polyhedron whose K 
pairs of faces are glued mathematically by holonomies. The defi- 
nition of fundamental cell is not unique, and in this paper by fun- 
damental cell we mean the Dirichlet domain seen from the "cen- 
ter" of the universe, D{X{)) - {x e [/ 1 Vy € F, |ji:-Xo| < Ijc-yjcol). 
In some space that is globally inhomogeneous, we can define 
a center, and the fundamental cells for 17 multiconnected flat 
spaces are found in Figure 1. No two points in a Dirichlet do- 
main can be linked by any holonomies; in this sense a Dirichlet 
domain represents the whole universe. In a globally homoge- 
neous space, the shape of the Dirichlet domain is independent 
of the observer's position, since it is always the same as that of 
the fundamental cell. In a globally inhomogeneous space, on the 
other hand, in general we do not stand at the center xq, so the ob- 
served Dirichlet domain D{Xobs) can vary from the fundamental 
cell D{xq). 

Because of the repeated structures given by holonomies, we 
can observe multiple copies (ghost images) of single objects in 
a multiconnected space that is small compared to the observed 
region, and most methods for constraining cosmic topology are 
based on this prediction. If we find such copies in the sky, it 
suggests that our Universe is multiconnected, while finding no 
ghosts indicates that our Universe is simply connected, at least 
within our sight, so the lower limit to the size of our Universe is 
obtained. 

There are mainly two methods for constraining cosmic topol- 
ogy. One is to use the 3D distributions of astronomical objects 
such as galaxies, quasars, and galaxy clusters, often called cos- 
mic crystallography method, treated in detail in this paper Use 
of the terminology of "cosmic crystallography" varies among 
authors, and in this paper we use it as the generic name for the 
3D data methods. This method can be classified into two types 
according to the type of topological signatures of spatial pat- 
terns they search for Those which search for type-I pairs (or, 
in general, n-tuplets), {x,y) and {yx,yy), with the relation of 
\x — y\ - \yx - yy\, have been proposed by, e.g. Roukema (1996) 
and Uzan et al. (1999). Roukema (1996) has developed a method 
of directly searching for two «-tuplets whose spatial distribu- 
tions are the same within a tolerance, while Uzan et al. (1999) 
have developed a method to count up the number of same sepa- 
rations seen in a catalog, which is expected to be larger in a small 
universe. Type-I pairs are found in every multiconnected space 
since any holonomies are isometrics that preserve distance. The 
others, which search for type-II pairs, {x, yx) and {y,yy) with the 
relation of \x — yx\ = \y — yy\, have been proposed by Lehoucq 
et al. (1996) and Marecki et al. (2005). Type-II pairs are only 
found in spaces where some of their holonomies are Cliff'ord 
translations, which translate all the points by the same distance. 
Lehoucq et al. (1996) have developed a method to detect Clifford 
translations using pair-separation histograms. The existence of 
such holonomies can be found as sharp spikes in histograms 
at the separations corresponding to their translating distances. 
Marecki et al. (2005) have improved the method by also using 
vectorial condition, not just separation. This "type-I and type-II 
pair" terminology was first introduced in Lehoucq et al. (1999). 

The other is to use the 2D cosmic microwave background 
(CMB) maps. Starobinsky (1993), and Stevens et al. (1993) have 
simulated CMB maps and angular power spectra of tempera- 
ture fluctuations, assuming 3-torus topology, and compared them 
with the observation by COBE satellite to obtain lower limits 
of the size of our Universe. Such simulations have been car- 
ried out by various authors (e.g. Luminet et al. 2003; Uzan et 
al. 2004; Riazuelo et al. 2004; Aurich et al. 2008; Aurich & 
Lustig 2010). These simulations and the observed quadrupole 



suppression suggest that we live in a small universe. There is 
another, direct observational method to detect topological sig- 
natures in CMB maps called circles-in-the-sky method (Cornish 
et al. 1998). This method uses pairs of circles with the same 
temperature fluctuation pattern, which may be the intersections 
of the last-scattering surface (LSS) and the observer's Dirichlet 
domain. Cosmic crystallography method lost its popularity af- 
ter the WMAP data release, since these CMB-based methods 
can constrain more topologies than the former, simply because 
the CMB data covers the larger region. Several authors have 
searched for these circles and obtained diverse results. Some 
claim that most of nontrivial topologies are ruled out (Cornish 
et al. 2004; Key et al. 2007). They have searched for antipo- 
dal or nearly antipodal pairs of circles in the WMAP map, and 
found no such circles to obtain the lower limit of the cell size 
as ~ 24 Gpc. This constraint cannot be applied to those spaces 
whose matched circles can be highly deviated from antipodal. 
Other authors claim that they have found the hints of multicon- 
nected spaces using improved versions of the circles-in-the-sky 
method (e.g. Roukema et al. 2008; Aurich 2008), which is incon- 
sistent with the former claim. This disagreement suggests there 
are methodological problems, which motivates us to revisit the 
cosmic crystallography method. 

As mentioned above, cosmic crystallography method has 
many versions since the original ones (Lehoucq et al. 1996; 
Roukema 1996) were introduced. These existing versions, how- 
ever, are no longer useful in universes that are comparable to the 
observed region in size due to the lack of the topological ghosts, 
which makes them unrealistic for practical application. This in- 
dicates the need to construct a more sensitive method, and in this 
paper we present such a new method that assumes Euclidean ge- 
ometry. We acknowledge that we have abandoned some gener- 
alities so as to enhance signal, but it is not serious because we 
have full mathematical knowledge about flat topologies. This pa- 
per focuses on methodology, and the simulated catalogs of "toy 
quasars" used here ignore some actual physics, such as spatial 
correlations and cycles of activity, since these effects do not af- 
fect the general results. 

In section 2 we review the 18 3D spaces described by 
Euclidean geometry and their holonomies. In section 3 we de- 
scribe the basics of the cosmic crystallography method and our 
new techniques for contrasting topological signals. To see the 
success of our improvements, we apply our method to catalogs 
of toy quasars. In section 4 we describe our catalogs and models 
of topologies to be tested. In section 5 we show the results and 
give some discussions. Conclusions are given in section 6. All 
calculations throughout this paper are done in comoving space. 



2. Flat spaces and their holonomies 

2.1. Classification of tlie 3D flat manifolds 

Throughout this paper we are interested in spaces with zero cur- 
vature, such as those, described by Euclidean geometry, as sug- 
gested by many observations (e.g. Qtot = 1 .0050;^QQggj from 
WMAP-hBAO+SN data, by Hinshaw et al. 2009). We therefore 
review a complete set of the 18 flat manifolds with three di- 
mensions. There are no other such topologies mathematically 
(Nowacki 1934), meaning that our Universe is one of them, if 
Qtot is exactly equal to unity. The classification is summarized in 
Table □ 

According to this table, most of the spaces are globally in- 
homogeneous. The shapes of their Dirichlet domains depend on 
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Symbol 


Name 


Compact 


Orientable 


Homogeneous 


El 


3-torus 


yes 


yes 


yes 


El 


half-turn space 


yes 


yes 


no 


Ei 


quarter-turn space 


yes 


yes 


no 


E4 


third-tuiTi space 


yes 


yes 


no 


£5 


sixth-turn space 


yes 


yes 


no 


Ee 


Hantzsche-Wendt space 


yes 


yes 


no 


El 


Klein space 


yes 


no 


no 


Ei 


Klein space with horizontal flip 


yes 


no 


no 


E9 


Klein space with vertical flip 


yes 


no 


no 


E]Q 


Klein space with half turn 


yes 


no 


no 


E\i 


chimney space 


no 


yes 


yes 


En 


chimney space with half turn 


no 


yes 


no 


En 


chimney space with vertical flip 


no 


no 


no 


£14 


chimney space with horizontal flip 


no 


no 


no 


El, 


chimney space with half turn and flip 


no 


no 


no 


Eli 


slab space 


no 


yes 


yes 


En 


slab space with flip 


no 


no 


no 


Eli 


simply connected space E'' 


no 


yes 


yes 



the observer's location, and especially their linked faces can be 
highly deviated from antipodal. Moreover, E(„ Ej, E^, Eg, and 
£10 have nonantipodal faces linked by holonomies that are inde- 
pendent of the observer's location in the universe (e.g. Riazuelo 
et al. 2004; Mota et al. 2010). Negative results to a multicon- 
nected Universe obtained by Cornish et al. (2004) and Key et 
al. (2007) using WMAP data are based on the searches for an- 
tipodal or nearly antipodal pairs of circles, so their constraints 
L > 24 Gpc cannot be applied to these spaces anyway. It means 
that these spaces remain as possible models of our Universe, and 
they may be detectable by a cosmic crystallography method. 

The 17 manifolds £1 , ■ ■ ■ , £17 are all multiconnected, whose 
fundamental cells are illustrated in Figure [1] The doors or other 
marks on the faces represent the means for topological gluing. 
The unmarked faces are glued straight across by parallel trans- 
lation, except for the Hantzsche-Wendt space, where six pairs of 
faces are glued in the same way as the doored one. 



2.2. Elemental transformations in flat spaces 

Hereafter, we use a 4D coordinate system (w, x, y, z) to consider 
the 3-plane w = 1 as the simply connected Euclidean space, so 
a point (x,y,z) in E-' is represented as {l,x,y,z). With this no- 
tation, transformations in 3 -Euclidean geometry can be written 
as 4 X 4 matrices. We introduce nine elemental transformations 
below. They are not commutative, but generators of the holon- 
omy group of each space can be derived by their products. In 
this sense, they are elemental. 

1. Parallel translation 




SL\lh Turn Spaw; "ITiird Turn S|iat;e Hanlzsche-Wendl SptKe wiih Maiical Flip with Half Tlim 




ChiimeySiiacewith Chiraiey S|nce *ith Cbimi*^ Si«ce wiUi 

\trticidHip Horizaital l-lip Half Tmiandllip 



Fig. 1. Illustration of fundamental cells for the 17 multicon- 
nected spaces with Euclidean geometry. The chimney space fam- 
ily has two compact dimensions, while the slab space family 
only has one compact dimension. The others have three compact 
dimensions; they are compact. The simply connected space E^* 
is not shown here. Courtesy of Jeffrey Weeks and Adam Weeks 
Marano. This figure was first published in Cipra (2002). 
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1 



These are the generators of parallel translation along three 
coordinate axes, by a distance L. Parallel translation induces no 
global inhomogeneity because it operates all points equally. 



TAE) ^ 



C 1 








1 





1 














1 





I L 








1 J 



2. Rotation 



OAB) = 



1 














1 














COS 6 


- sin^ 








sin 9 


cos 6 



0,(9) = 



1 














COS 





sin 








1 








- sin 





cos 6 



3 



H. Fujii & Y. Yoshii: An improved crystallography method for flat spaces 



OM = 



1 














COS 9 


— sin 9 








sinfl 


cos 9 














1 



These are the generators of rotation about three coordinate 
axes, by an angle 6. Rotation is a globally inhomogeneous isom- 
etry. For example, a point near the x-axis is translated to a 
nearby point, while a point distant from the x-axis is translated 
to a distant point, by the same rotation Ox(ff). 

3. Reflection 



] 






C 1 




















1 J 






I 


( 1 








^ 





1 














1 





I 








-1 J 



These are the generators of reflection about three planes per- 
pendicular to coordinate axes, e.g. is a reflection about the y-z 
plane. Reflection is a globally inhomogeneous isometry. For ex- 
ample, a point near the y-z plane is translated to a nearby point, 
while a point distant from the y-z plane is translated to a distant 
point, by the same reflection R^. Moreover, the nonorientability 
of a nonorientable space comes from reflection. 

2.3. Neighboring hoionomies 

With these elemental transformations, generators of the holon- 
omy groups of each space are given as their multiplications. 
Specifically, neighboring hoionomies, i.e., hoionomies that link 
an object in the observer's Dirichlet domain to its topological 
copies in neighboring cells, are also written as products of ele- 
mental transformations as shown in Table |2] In sections 4 and 5 
we consider the most conservative cases in which neighboring 
hoionomies are barely detectable. 

Along with Table 1, it can be seen that the spaces whose 
holonomy groups contain reflection or rotation are globally in- 
homogeneous, and those whose holonomy groups contain reflec- 
tion are not orientable. 

Neighboring hoionomies can be classified into six types; 
half-turn corkscrew motion, quarter-turn corkscrew motion, 
third-turn corkscrew motion, sixth-turn corkscrew motion, and 
glide reflection. Hoionomies that come from different topolo- 
gies, but belong to the same type, are indistinguishable by them- 
selves. For example, if we observe a pair of antipodal transla- 
tions, we live in one of the spaces that have translation as neigh- 
boring ones, but we cannot specify in which space we live. 

In the following, the spaces in the parentheses are those that 
have the type of hoionomies as neighboring ones. 

1. Translation 

(E\,E2,E3,E4,Es,Ej,Eii,Eg,Eio,Eii,Ei2,Ei3,Ei4,Ei(,) 
Translation induces antipodal faces of the Dirichlet domain 
that are glued straight across. Translation is subdivided into 
three types of one-pair translation (type I), two-pair translation 
(type II), and three-pair translation (type III). Two-pair trans- 
lation only exists in quarter-turn space Ey, its four hoionomies 
Txi+L) and Ty(+L) have the common translating distance L, so 
either all four hoionomies are detectable or all of them are un- 
detectable. Three-pair translation, on the other hand, exists in 



third-turn space E4 and sixth-turn space £5; their six hoionomies 

T,{+\L)Ty(±^L) and T,(+L) have the common translating dis- 
tance L, so either all six hoionomies are detectable or all of them 
are undetectable. One-pair translation exists in the other spaces. 

2. Half-turn corkscrew motion {E2, E(,,Eia, £12, -Eis) 
Half-turn corkscrew motion is the half-turn followed by par- 
allel translation. Half-turn corkscrew motion is subdivided into 
two types of those that have rotational axes corresponding to 
translational directions (type I), and those with different rota- 
tional axes and translational directions (type II). Only E(, has the 
latter type, while the other four spaces have the former 

3. Quarter-turn corkscrew motion (£3) 

Quarter- turn corkscrew motion is the quarter- turn followed 
by parallel translation. Only £3 has quarter-turn corkscrew mo- 
tion, so detecting such hoionomies suggests that we live in a 
quarter- turn space. 

4. Third-turn corkscrew motion (£4) 

Third-turn corkscrew motion is the third-turn followed by 
parallel translation. Only £4 has third-turn corkscrew motion, so 
detecting such hoionomies suggests that we live in a third-turn 
space. 

5. Sixth-turn corkscrew motion (£5) 

Sixth-turn corkscrew motion is the sixth-turn followed by 
parallel translation. Only £g has sixth-turn corkscrew motion, so 
detecting such hoionomies suggests that we live in a sixth-turn 
space. 

6. Glide reflection (£7, £», £9, £10, £13, £14, £15, £17) 

Glide reflection is the reflection followed by parallel transla- 
tion. They are subdivided into two types of those that have reflec- 
tional planes parallel to translational directions (type I) and those 
with nonparallel reflectional planes and translational directions 
(type II). The six spaces of £g, £9, £0, £14, £ig, and £17 have 
the former type, while the four spaces £7, £8, £9, and £10 have 
the latter. 



3. Method 

3.1. Basic ideas 

Cosmic crystallography is a series of statistical techniques to ex- 
tract topological information from a given astronomical catalog 
first proposed by Lehoucq et al. (1996) and Roukema (1996). 
The former method is to search for type-II pairs, and the latter 
is to search for type-I n-tuplets (generalization of type-I pairs to 
n objects). This type-I, type-II terminology is that of Lehoucq 
et al. (1999), and is unrelated to the definitions of types in 
Section 2.3. In this paper we propose a new version of crystal- 
lographic method, which is mainly for collecting type-I pairs as 
in Roukema (1996) and Uzan et al. (1999) and has filters that 
are related to those of Marecki et al. (2005). Flat spaces are as- 
sumed throughout the paper, but the methodology in this section 
can be applied to spaces with any curvature. 

We assume that our Universe has topology M = /F, where 
is a simply connected 3-Euclidean space, and F is a holon- 
omy group on E^*. Suppose that in our catalog we have objects 
Pi, ■ ■ ■ ,Pn, whose comoving positions are given by Xi, ■ ■ ■ ,Xn, 
respectively. We often indicate the object P,- itself by its posi- 
tonal vector jc, for convenience. Ghosts of two objects Xj and Xj 
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Table 2. Neighboring holonomies of each topology represented as products of elemental transformations. 



Space symbol 


Neighboring holonomies 


El 


lA±Ei), Jy{±E2), Jz(±E3) 


El 


y.t(±Li), ly(±L2), l,(±L3)0-{n) 


E3 


Jx(±Ei), ly(±Li), 1 ,(±L2)0-(±n/2) 


£4 


TMLi)Ty(± V3Li), T,(±2Li), T,(±L2)0-A+2n/3) 


£5 


TA±Li)Ty{± V3Li), r,(±2Li), T,(±L2)0,(±n/3) 


E6 


TA±Li)T,{±L2)0An), Ty(±L2)T,{±L3)0,(7x), r,(±L3)r,(±L,)0,(;r) 


El 


TA±Li)Ty{±L2)Ry, r,(±2Li), T,{±Li) 


Es 


TA±Li)Ty{±L2)R„ TA±2Li), U±L3)R, 


E9 


TA±Li)TA±L2)Ry, TA±2Li), T-S±L3)R,. 


E\o 


TA±Li)Ty{±L2)Ry, r,(±2Li), TA±Li)0,{7:) 


Ell 


lA±Ei), ly{±L2) 


En 


TA±Li), T,(±L2)0,{n) 


En 


TA±Li), T,{±L2}R, 


El4 


TA±Li), T,(±L2]R, 


En 


TA±Li)R,, r,(±L2)0,(;r) 


E16 


U±Li) 


En 


T-A±Li)R, 


En 


none 



by a holonomy y e F are given by yjc, and yxj, respectively. All 
holonomies are isometrics that preserve the distance, so we have 

\Xi-Xj\^\yXi-yXjl (2) 

independently of the holonomy's other properties (Figure Q. If 




Fig. 2. Two pairs of objects linked by a holonomy y ideally have 
the same separations. This ideal situation breaks down owing to 
peculiar velocities and observational limits. 

we find two pairs {Xi,Xj) and (jc^., xi) such that jCjl - \xk-xi\, 
they can be ghosts of each other. We call the two pairs of objects 
a quadruplet hereafter (notice the difference between the type- 
I n-tuplet and the quadruplet here), and the statement above is 
equal to saying that a quadruplet [(Xi,Xj),{Xk,Xi)] with the re- 
lation \Xi - X i\ = \xii - xi\ can be topological, and finding such 
quadruplets indicates the multiconnectedness of our Universe. 
However, in practice, these separations cannot be exactly the 
same because of various effects such as peculiar velocities and 
observational limits, so we search for pairs having separations 
deviated at most by some positive value s. Because of this finite 
spatial resolution, some nontopological quadruplets stochasti- 
cally pretend to be topological. The number of such quadruplets 



M is expected to be proportional to the total number of quadru- 
plets X 4C2 X ^, so the topological index is given by 

S - ^ - 8^ (3) 

^C4X4C2X^ A^(A^-l)(A^-2)(A?-3)' 

In a small universe, there are real topological quadruplets added 
to false stochastic ones, so the index S will take a higher value 
than expected in a simply connected universe. This formulation 
is slightly different from the CCP method, but their essence is 
common to both. 

Until now we have only used the general property of 
holonomies as isometrics, i.e. preserving the distance. If a holon- 
omy y is a Clifford translation, which translates all points by the 
same distance, additional pairs have the same separations: 

\Xi -yXi\^\Xj-yxj\. (4) 

In this case, the topological signal by y is multiplied by a factor 
of 2. Since we concentrate on spaces with Euclidean geometry, 
a Clifford translation is a normal parallel translation. 

3.2. Sophistication 

The existing crystallography method has been verified as use- 
ful for constraining cosmic topology, only for topologies whose 
fundamental cells are so small that the observed region contains 
no fewer than several tens of cells within it. However, in reality, 
the fundamental cell of our Universe can be comparable to the 
observed region in size. Of course it can be larger than the ob- 
served region, but in this case no methods can be used, so we do 
not consider it. For such situations, the crystallography method 
is no longer valid because the real signal is too weak compared 
to the false signal, which occurs by chance. It is necessary for us 
to make the methods more sophisticated in order to decrease the 
false signal, while keeping the real signal. 

3.2.1 . Imitation of tine circles-in-tlie-sky metliod 

First of all, we propose to concentrate on a thin, shell-like region 
ri < r < r2, like the circles-in-the-sky method using LSS. The 
reason we consider such a region is that we have to reasonably 
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decrease the number of objects, since, roughly speaking, the real 
signal is proportional to A^^, while the false signal is proportional 
to A^'*. If we can detect the signal of a space that is comparable 
to the shell region in size, then any smaller spaces are detectable 
since they produce stronger signals. In principle, therefore, the 
larger the region we consider, the more topologies we can con- 
strain. However, the number of objects becomes too small when 
the region is too large, so we have to find a good balance. 

When the region is chosen, then we have to put all the 
quadruplets through "filters" to drop any false stochastic ones. 
Here we introduce three filters: the first one has already been 
discussed both in the previous section and in Uzan et al. (1999), 
but the others are new ones. 

3.2.2. The first filter : separation 

All holonomies are isometrics that preserve distance, and we 
consider the condition 

\\Xi-Xj\-\Xk-Xi\\<E (5) 

as the first filter to extract information about topological tensing. 
This corresponds to selecting type-I pairs as in Roukema (1996) 
and Uzan et al. (1999). This filter, because of the finite resolu- 
tion e, is insufficient when the real signal is weak. Many false 
quadruplets pretending to be topological hide the real ones. This 
filter can be used in spaces with any curvature, unlike the other 
two, since it uses the general property of holonomies. 

3.2.3. Tine second filter : vectorial condition 

To extract information about repeated structures due to multi- 
connectedness, we should pay attention not only to separation 
but also to vectorial properties of holonomies. Vectorial prop- 
erties are not common among holonomies, unlike the distance- 
preserving property, so this filter cannot constrain all holonomies 
at once. Here we introduce five kinds of filters that constrain their 
specialized holonomies. These choices are not unique; in partic- 
ular, we can construct more specialized filters to constrain more 
specific holonomies, but it induces a drastic increase in free pa- 
rameters. In some sense, our five filters can be regarded as the 
generalizations of bunch-of-pair (BoP) selection in Marecki et 
al. (2005), which selects type-II translational pairs that have both 
the same separations and the same directions, to nontranslational 
holonomies. 

1 : Translation 

This filter is for constraining the parallel translation y - 
T^^{Li)Ty{L2)T~(L3). When a given quadruplet [{Xi,Xj),(Xk,Xi)] 
is a y-quadruplet, there are two possibilities such that yixi, Xj) - 
(Xk, Xi) or yiXi, xj) - {xi,Xk)- In the former case, we ideally have 
a relation: 

Xi - Xj ^ Xk- xi, ji - jj ^jk- yi, Zi - Zj ^ Zk- zi- (6) 
In the latter case, on the other hand, we ideally have a relation: 
Xi-Xj = -{xk-xi), yi-yj = -(yk-yi), Zi-zj = -(zk-zi). (7) 
We consider these possibilities together to construct a filter: 

\(xi- Xj)±(xk- xi)\ < (8) 

\(yi-yj) + {yk-yi)\<Sy, (9) 
\(zi-Zj)±(zk-z,)\<e„ (10) 



independent of the translating distance and direction. Note that 
s = (e^ , , e ) here is different from s in the first filter This fil- 
ter corresponds to the BoP selection in Marecki et al. (2005), 
but our filter also detects type-I pairs, not only type-II pairs. 
However, this means that our filter also detects false type-I and 
type-II pairs, so there seems to be no significant advantages and 
disadvantages between them. 

2 .• Half-turn corkscrew motion and glide reflection 

This filter is for constraining half-turn corkscrew motion, e.g. 
yi - Tx{L\)Ty{L2) T-{L2)0-{7t), and glide reflection, e.g. 72 = 
T^{Li)Ty{L2)T,{L3)R-. For a given quadmplet [(Xi,Xj),{Xk,Xi)] 
to be a 7„ -quadruplet, there are two possibilities such that 
yniXi,Xj) = ixk,xi) ory„iXi,Xj) = ixi,Xk). 
In the former case, we ideally have a relation for n - I 

Xi - Xj = -(xk - xi), yi -yj ^ -(yk - yi), Zi - zj ^ Zk- zi, (11) 
and for « = 2 

Xi - Xj ^ Xk- xi, yi - yj = yk - yi, Zi - Zj = -(zk - Zi). (12) 

In the latter case, on the other hand, we ideally have a relation 
for n - I 

Xi - Xj ^ Xk- xi, yi - yj ^yk- yi, Zi - Zj = -(Zk - Zi), (13) 
and for « = 2 

Xi -Xj ^ -{xk - xi), yi -yj ^ -(yk - yd, Zi - zj ^ Zk- zi- (14) 

We consider these possibilities together to construct a filter: 

\(xi- Xj)±(xk- xi)\ < s^, (15) 

\(yi-yj) + (yk-yi)\<Sy, (16) 

\(zi-Zj) + (zk-zi)\<s,. (17) 

Here the rotation and reflection are based on the "global" di- 
rections of the universe about which we primarily have no idea, 
so in practice we have to find them by changing our choices of 
coordinate axes (see section 5.2.1 for details). 

3 : n-th turn corkscrew motion (n —4,3,6) 

These filters are for constraining n-th turn corkscrew mo- 
tion for n - 4,3, and 6. It is, for example, a holonomy 7 = 
T^{Li)Ty{L2)T.{Li)0~{2n/n). As in the previous filters, there are 
two ways to pass through this filter for a given quadruplet. 

We construct the filters by considering the two possibilities 
together, say, 

\{xi - Xj) cos{2n/n) - (yi - y j) sin(27r/n) ± (xk - x/)| < s^, (18) 
\(xi - Xj) sin{2n/n) + (yi - y j) cos(27r/n) + (yk- yi)\ < Sy, (19) 
\(Zi - Zj) + (zk - zi)\ < s,. (20) 

Corkscrew motion includes rotation based on the global direc- 
tions, so we have to change our choices of coordinate axes to 
find it, as already mentioned. 

An additional notice here is that in our calculations 
we do not distinguish two quadruplets [(Xi,Xj),(Xk,Xi)] 
and [(Xk,Xi),(Xi,x j)], so the above filters will ignore 7"'- 
quadruplets. For this we also use the filters 

\(xi - Xj) cos(27r/«) + (yi - yj) sin(2n/n) + (xk - xi)\ < Sf,(2l) 
I - (xi - Xj) sin{2n/n) + (yi - y j) cos(27r/«) + (yk- yi)\ < £,(,22) 
\(z,-Zj)±(zk-z,)\<s„ (23) 

which are unnecessary for the previous two cases since 7"'- 
quadruplets can also be constrained by the same filters. 
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3.2.4. The third filter: lifetime of objects 

Astronomical objects have finite Ufetimes fiifg. Suppose two ob- 
jects Xi and Xj whose cosmic times are f, and tj, respectively. In 
order that they are ghost images of each other, 



< hife 



(24) 



is a necessary condition, so a quadruplet [(Xi,Xj),(Xk,Xi)] should 
be dropped unless Af;^, Afy/ < fiifg or Atii,Atjk < fiife. This very 
closely corresponds to the redshift filter in Marecki et al. (2005), 
and the difference is only an expression. Some preceding studies 
ignore this effect and use all pairs despite their cosmic times. 
This is a crucial fault unless L <sc cfufe, where L is a characteristic 
size of the universe. 



3.2.5. Finishing: classification of objects 

Though the above three filters can drop false stochastic quadru- 
plets while keeping real topological ones, it is possible that the 
false signal is still strong enough to hide the real signal. For this 
situation, we introduce an additional technique to contrast real 
signals. 

Each object Xj [i = I, - ■ ■ ,N) is assigned an integer i,, the 
number of final candidate quadruplets that have passed through 
all filters and include jc, as their members. False signals are 
stochastic, so the members of false quadruplets are randomly 
distributed. As a result, the possibility that many final candidates 
share a single, common object is small, so i, rarely takes a high 
value. Real signals, on the other hand, all come from ghost pairs. 
If there are n y-pairs (jc,-, yXj), ■ ■ ■ , (xj, yx j) where y is not a par- 
allel translation, > n - 1 for each x^, whereas > 2(« - 1) 
if y is a parallel translation. As a result the histogram of 5, 
will contrast real topological signals from false stochastic ones. 
This technique corresponds to performing the «-tuplet search- 
ing in Roukema et al. (1996), for all n simultaneously. The CCP 
method by Uzan et al. (1999) has no such techniques, and is 
therefore unable to distinguish, for example, one real «-tuplet 
from n(n - l)/2 false pairs. 

Through these techniques, candidates of topological ghosts 
are sampled in a catalog, and their distributions in the sky give us 
a hint of the topology of our Universe. They trace the topological 
gluing of the faces of our Dirichlet domain. 



3.3. Spatial resolutions 

In this paper we do not consider any technical uncertainties; i.e., 
we assume that we have the full-sky, full-redshift, and 100% 
complete catalog with infinite accuracy. However, even in such 
a situation, peculiar velocity with respect to the comoving flame 
prevents us from having the ideal catalog where Eq.(2) holds 
precisely. Effects of peculiar velocities consist of two parts: (i) 
integrated effect and (ii) instantaneous effect. They were first dis- 
cussed by Roukema (1996), while these terminologies are given 
by Lehoucq et al. (2000). The integrated effect comes from the 
true motion of objects, and the instantaneous effect comes from 
the fake positional change due to an additional redshift by the 
Doppler effect. To see these effects quantitatively, we consider 
two objects Pi and Pj = yP,, linked by a holonomy y. Their 
positions are given by x, and Xj - yXi, respectively, when the 
object is comoving. In practice, however, the observed positions 
are slightly different from those. 



First we consider the integrated effect. Denoting the cosmic 
times of P, and Pj by f, and tj, respectively, we obtain 



r'j V 
' J,, a(t) 



(25) 



as a motion of the object from f = f,- to f = tj, so it affects 
the observed position of Pj as an additional term to yjc,. Here 
Aint^:, = by definition. 

Next we consider the instantaneous effect. Assuming the pe- 
culiar velocity v = |v| is time-independent, the extra redshift by 
the Doppler effect is given by 



AZdop = —(1 +Zcos), 

c 



(26) 



to first order in vy/c, where v\\ is a radial component of v, and 
Zcos is a purely cosmological redshift. This extra term affects the 
observed comoving radial distance of objects via 



Ains'" 



^0 



da 



(27) 

This one, unlike the integrated effect, affects both P, and Pj on 
their observed positions. 

Incorporating these two effects simultaneously, the observed 
positions of P, and Pj are given by 



r°''" -X- + A- r — 

\Xi\ 



(28) 



and 



Xf = yXi + AintXj + AinsO— — ;— (29) 



\yxi + AintXjl 



respectively. Here we have seen the effects of peculiar veloci- 
ties on observed positions by which we would fail to recognize 
real topological quadruplets, which helps us to determine spatial 
resolutions s in the first filter and e in the second filter. They 
should be low enough to avoid missing the topological quadru- 
plets that are distorted by these effects, however, simultaneously, 
they should be high enough to decrease the false signal. A bal- 
ance between them determines s and s (details are in section 5). 



4. Simulations 

4. 1 . Toy catalogs of quasars 

To test our new method, on the level of methodology, we gen- 
erated toy quasar distributions in flat spaces with the stan- 
dard A-CDM cosmology (n,„ = 0.27, Oa = 0.73, i/o = 71 
km/sec/Mpc), and applied the method to them. We make com- 
plete, full-sky and full-redshift catalogs with infinite accuracy 
(more precisely, double-precision floating-point accuracy). As 
the nature of quasars, we consider the cosmological evolution of 
comoving density, peculiar velocity, and lifetime, since the most 
crucial factor is a ratio of the number of topological quadruplets 
to that of the false quadruplets. The lifetime of quasars concerns 
the former, and the comoving density and peculiar velocity con- 
cern the latter The other natures of clustering and cycles of ac- 
tivity are ignored because we focus on methodology here, and 
these effects are not important at this level. The method itself 
can be applied to catalogs of any extragalactic objects, such as 
active galactic nuclei (AGNs), high-z galaxies, and galaxy clus- 
ters. 
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The redshift evolution of the comoving density of quasars is 
taken from Osmer (2004), but we have extrapolated it to z ~ 7. 
We simplified the luminosity evolution of quasars such that they 
emit radiation with constant luminosity during the fixed duration 
hife -10^ yr. We also simplified the peculiar motion of quasars 
to move with constant speed v = 500 km/sec and with randomly 
chosen directions. 

Every 10^ yr from f = to the present f = 1 .37 x 10'" yr, new 
quasars are generated and randomly distributed, and the quasars 
reaching their lifetime of 10^ yr are removed. For simplicity, no 
spatial correlations are considered in this paper While the co- 
moving density of quasars used in simulations is common, the 
actual number of quasars listed in the catalogs is different. This 
can be viewed as the cosmic variance. 

The topological signal is searched for in a shell region 7.8 
Gpc< r <8.2 Gpc (corresponding to 4.7 < z < 5.5), in which 
about 2000 quasars are contained. This choice of radial width 
does not necessarily have significant meanings, but this is close 
to the limiting size {z ~ 7) of the actual distribution of quasars. 
Hereafter, any catalog mentioned is referred to as a part of the 
full catalog satisfying 7.8 Gpc< r <8.2 Gpc. 

4.2. Models of spaces 

For each topology, the holonomies we target are assumed to have 
translational parts, e.g., T,{Lt,) for T,{L^)0-{ji), satisfying 

+Ll + Ll = L = 16 Gpc. (30) 

which is a limit size that can be detected with our catalog 7.8 
Gpc< r <8.2 Gpc. If they are detectable, then any other space 
that is smaller than them is also detectable. It would be difficult 
to make a detection in this case using other 3D methods. 

In each calculation, our policy is "holonomy first". 
According to the policy, unless otherwise stated, in each case 
we consider the most conservative cases in which all holonomies 
except for the targeted ones are beyond the observed region: they 
were undetectable. In other words, we consider the severest case 
for each topology to detect it. For example, when we target the 
z-parallel translations T^(±L) of 3-torus, we neglect the other 
translation Tx{+L\) and Ty(+L2) by considering L\,L2 oo. 
When we target the parallel translation of third-turn space, on 
the other hand, we must simultaneously consider the translation 
in all directions because of their common translating distance, 
(see Table 2) 

5. Results and discussions 

5.1. Improvements in detecting topological signal 

Here we show the results in detail for one type of holonomy. 
Those for the other types are given in the appendix. Two catalogs 
of toy quasars, one in a slab space with Tr{+L) and the other in 
the simply connected Euclidean space with no holonomies, are 
compared. Each catalog contains 1980 quasars. Among the 1980 
ones in the former catalog, there are six pairs of ghosts, meaning 
that there are 30 (=6^2 x 2) topological quadruplets. They are 
hidden in the large number of false stochastic ones. To see the 
effects of three filters introduced in section 3, we have increased 
the number of filters one by one, and counted the number of 
quadruplets that have passed through the filters. 

As for spatial resolution, we used e = 3.7 Mpc in the first 
filter and e - (6.0, 6.0, 6.0) in the second filter in units of Mpc. 
They are determined as follows. We generated 10000 topological 



quadruplets [(Xi,Xj),{T^(+L)Xi,T^(+L)Xj)]. The value of s was 
chosen such that 99% of 10000 quadruplets satisfied the condi- 
tion -jCjl -|r;(+L)j(:,- r-.(+L)j!:y|| < £. Then, each component 
e,. By, and were determined in the same way as e. We denote 
s' = maxlsv, Sy, £;) and conservatively substitute {s',e',s') for 
s. In principle, we have to set s and s such that 100% of the 
quadruplets satisfy the condition, but it needs very high resolu- 
tion, and meaninglessly increases the false signal. We balanced 
them to obtain the values. 

The total number of quadruplets is given by 

M = 1980C4 X 4C2 X ^ = 1915375615065. (31) 

This drastically large number of quadruplets are filtered, and the 
false stochastic ones are dropped as shown in Table |3] We can 
find that the first filter (separation) is insufficient for detecting 
holonomies when there are only a few ghosts in this case. It 
suggests that the existing crystallography method like the CCP 
method is no longer useful for this situation, since the real signal 
is too weak. The second filter plays an important role in extract- 
ing the real signal. The third filter only plays a compensative role 
here, but will be important in somewhat more delicate situations. 



Table 3. The number of quadruplets (A^quad) that have passed 
through the filters of (1) separation, (2) vectorial condition, and 
(3) lifetime of objects. 



Filtering condition 


A^quad in the 
slab space 


A^quad in the simply 
connected space 


no filters 


1915375615065 


1915375615065 


filter (1) 


663853290 


664887124 


filter (l)-i-(2) 


342 


323 


filter (1)-H(2)-H(3) 


253 


239 



Though the false signal has been decreased successfully, it 
still hides the real one, so we additionally need the classification 
of objects (section 3.2.5). The i,-histogram, where Sj is the num- 
ber of final candidate quadruplets, which include the object P,-, 
for each catalog that is fully filtered is shown in Figure |3] There 
is a hill that peaks at i, - 9, in the histogram for the slab space. 
It may be constituted by topological ghosts, as expected in sec- 
tion 3.2.5. To check this in Figure[3]we display the quasars with 
Si > 7 using the Lambert's azimuthal projection in the z-x plane. 
The pair of five objects are clearly linked by parallel translation, 
confirming that our space has a translational pair 

Other examples for this projection are given in Figure 
where we get an advance of the results in Appendix. Quasars 
with high Si values, which constitute the hills, are those of a one- 
pair translation, as in the previous figure, half-turn corkscrew 
motion (type I), and glide reflection (type I). They clearly trace 
the way for topological gluing for each type of holonomies. If 
the objects that constitute the hills are randomly distributed, it 
is likely that they are merely stochastic. Moreover, in this stage, 
even if there are still some false stochastic ones, it is possible 
to find repeated patterns of real signal with somewhat additional 
techniques, if exist. 

A relatively large number of quasars at i, - 2,4, and 6 in 
both histograms are due to a condition for quadrilaterals to be 
parallelograms, i.e., one pair of opposite sides that are parallel 
and equal in length. If a quadruplet [{Xi,Xj), {Xk,xi)] has passed 
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Fig. 3. Left: the i,-histograms of the two catalogs with 1980 
quasars in the simply connected Euclidean space and in the slab 
space. Si is the number of final candidate quadruplets that include 
the quasar P, . In the slab space, a hill constituted by topological 
ghosts is clearly seen. Right: The objects with > 7 in the slab 
space. It can be seen that they are linked by parallel translation. 



Fig. 4. The quasars that constitute the hill in the histogram for 
each simulation. Left: one-pair translation > 7), middle: half- 
turn corkscrew motion of type I (s, > 6), right: glide reflection 
of type I (i, > 6). They clearly trace the topological gluing for 
each holonomy. 



through all the filters, the quadrilateral PtPjPkPi is close to a par- 
allelogram, and then another quadruplet [{Xi,Xit),{Xj,Xi)] also 
tends to pass through the filters. This effect enhances the num- 
ber of quasars with even values of s/. Any other zigzag features 
of this and the other histograms (Figures l5T lA.2t are not mean- 
ingful, but are merely stochastic. 

5.2. Some complexities due to global inhomogeneity 

We have seen that our method is so sensitive to real signal that 
even the limiting case can be detected. A slab space topology has 
been assumed there, but the method is valid in any spaces as will 
be shown in Appendix. However, in globally inhomogeneous 
spaces (see Table 1), situations become somewhat complicated. 
First, globally inhomogeneous spaces have specific, preferred 
coordinate systems by which their holonomies are defined, and 
our filters are also based on them. It means that there is no signal 
unless we use the correct coordinate system. Second, the shape 
of the Dirichlet domain depends on the observer's location in the 
universe, so the real signal may not be seen for some observers. 
We study these effects below. 

To do this, we take, for example, a Hantzsche-Wendt space 



topology with Li = L3 



^/2 



11.3 Gpc, and L2 



since all the "abnormal" effects due to global inhomogeneity 
can be seen in this space. With these parameters, four half- 
turn corkscrew motions T~{+Li)T^{+Li)0.{7r) are should be de- 
tectable. 

As for spatial resolution, we have used e = 5.8 Mpc and 
e = (7.3, 7.3, 7.3) in units of Mpc. These values are determined 
in a similar way to that in the previous section, but we have set 
the lower confidence level of 90% and have taken the average 



values over 1000 observers randomly distributed in the funda- 
mental cell. With the confidence level 99%, spatial resolution 
becomes so high that the false signal dominates in some cases. 
Taking this into account, we proceed with the value of 90%. 

5.2.1 . Choices of coordinate axes 

In reality, we do not know the global directions of the Universe, 
so we have to find them, if they exist, by changing our choices 
of the coordinate axes. This situation is analogous to the circles- 
in-the-sky method where we have to search for the centers of 
circles. 

As long as the z-axis is correctly chosen, rotation about 
the z-axis does not affect anything, when we consider 
T-(+L/ V2)r,(+L/ y/2)0-{K), so we only have to find the z-axis 
with some accuracy. We can roughly estimate the angular reso- 
lution A6 by 

-LAO ~ s. (32) 
2 

Substituting L - \6 Gpc and e - 7.3 Mpc, we derive AO ^ 
0.05° ^ 3'. 

We have prepared two catalogs containing 1983 quasars in 
the Hantzsche-Wendt space in which the observer is at the center 
and in the simply connected space. The former catalog includes 
13 pairs of ghosts. Before searching for quadruplets and apply- 
ing our filters, we rotated the catalog about the y-axis from the 
initial state, by AO - V ,5' , 10', and 0.5°. The results are shown 
in Figure|5] It can be seen that the limit is in between 5' and 10', 
roughly corresponding to the pre-estimated value AO ^ 3'. The 
area of the half sky is Afuiisky/2 = 360^ /27r deg^, so the number 
of trials needed to find the z-axis is 



AMisky/2 _ 360V27r x 60^ 
(A0)2 ~ 



52 



3 X 10" 



(33) 



where we use AO - 5'. Simulation runs are mutually indepen- 
dent of each other, which takes about a few minutes with an or- 
dinary personal computer Therefore, a whole simulation can be 
accomplished faster by using many computers. No special com- 
puters are needed to clear up the grand mystery of the shape of 
our Universe. 



Euclidean space 
No rotations 
1 minute 
5 minutes 
10 minutes 
0.5 degrees 




Fig. 5. The i,- -histograms in the Hantzsche-Wendt space rotated 
about the y-axis by several angles. For comparison, the s,- 
histogram is shown in simply connected space. 
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5.2.2. Observer's location 

Next we consider cases in which the observer deviates from the 
center of the universe. We do not enter into the general theory 
of this effect, and simply show that the method works well even 
in these situations. We have located the observer at several posi- 
tions in the universe from which at least one pair of holonomies 
is detectable. 

For an observer at Z = {X, Y, Z), the condition \X - yX\l2 < 
r2 should be satisfied in order to detect a holonomy y in a shell 
region of ri < r < ri- The lefthand side is the distance from the 
observer to the face of his Dirichlet domain, which corresponds 
to y. Continuing with the same Hantzsche-Wendt space as an 
example, the lefthand side can be written explicitly as 




for y = r,(+L/ y/2)T_,{L/ ^l2)0-{n), and 




for y = T-i±LI y/2)TA-L/ y/2)0,(n), where L = 16 Gpc and 
r2 - 8.2 Gpc. An observer satisfying each condition can detect 
each pair of holonomies. It can be seen that the Z-location of 
the observer affects nothing. However, the bigger the F-location 
becomes, the farther away all faces are located, hence the harder 
to detect all holonomies. The X-location has an effect of pushing 
two faces away, while drawing the other two, which is due to the 
translational direction that does not accord with the rotational 
axis. The last effect cannot be seen in spaces without half-turn 
corkscrew motion (type 11) or glide reflection (type 11). 

According to these conditions, we chose six locations of 
{X, Y,Z) = (0, 0, 0), (1, 0, 0), (3, 0, 0), (5, 0, 0), (0, 0.5, 0), and 
(0, 1, 0), in units of Gpc. We prepared catalogs of toy quasars 
seen from these observers, and applied our method to them. The 
coordinate axes were chosen correctly here. Results for these ob- 
servers are given in Tableland Figure |6] 

In this table, it can be seen that the signal gets weaker as 
the observer moves along the y-axis, since the faces get farther 
and the ghosts decrease in number This is also seen in the his- 
tograms, where the hills constituted by ghosts disappear It is 
necessary to use a larger shell region, if possible, to detect the 
ghosts. 

As the observer moves along the x-axis, on the other hand, 
the signal and the hills remain like the initial one in which the 
observer is located at the center, since one pair of faces gets far- 
ther, but the others get closer As long as the coordinate axes 
are chosen correctly, our method can detect holonomies that are 
close enough to the observer Our method is suited to obtaining 
a lower limit to the size of the Universe. 



6. Summary and conclusions 

In this paper we have developed a cosmic crystallography 
method that significantly extends previous methods. A thin, 
shell-like part (ri < r < ro) of the full catalog is used, sim- 
ilar to the circles-in-the-sky method. This region should be as 
large as possible, but the comoving density of objects should 
be simultaneously high enough there. The quadruplets of these 
objects are filtered three times; separation, vectorial condition, 
and lifetime of objects. These filters drop the false stochastic 




Fig. 6. The dependence on the observer's location in the uni- 
verse. Top left:(0, 0, 0), top right:(0, 0.5, 0), mid left:(0, 1, 0), 
mid right:(l, 0, 0), down left:(3, 0, 0), down right:(5, 0, 0). 



ones, while keeping the real topological ones. Flat spaces de- 
scribed by Euclidean geometry are assumed, and the second fil- 
ter, vectorial condition, is specialized to detect their holonomies. 
This assumption is not fundamental, since similar filters for the 
holonomies in spherical or hyperbolic spaces can also be con- 
structed. 

The number of quadruplets having passed these filters is 
translated into an index of multiconnectedness of the universe, 
and then the objects are classified by s-,, the number of such 
quadruplets that include the object P, as one of their members. 
The 5, -histograms emphasize the topological ghosts: if ghost im- 
ages really exist, they show some hills in the large- i, region. 

We tested the method in catalogs of toy quasars with no ob- 
servational uncertainties. We considered the severest cases for 
each type of holonomies, and any holonomies except for the tar- 
geted ones are assumed to be beyond the observed region. The 
translational distance of each holonomy is fixed with L - \6 
Gpc for which 8 = (16/2) Gpc corresponds to z ~ 5, and the 
shell region is chosen as 7.8 Gpc< r <8.2 Gpc (4.7 < z < 5.5). 
In this situation, where the space is comparable to the observed 
region in size, existing crystallography method is no longer valid 
due to the contamination of false signal. It is found that our fil- 
ters are able to eliminate such contamination, and the existence 
of topological ghosts is clearly recognized by the i, -histograms. 
Our method considers isometries of flat spaces generically ac- 
cording to their type (translation, corkscrew motion, or glide re- 
flection), without requiring the specific manifold to be chosen 
a priori. This is loosely analogous to the generic nature of the 
circles-in-the-sky method. Except for the case of translation, we 
have to search over several million possible orientations of the 
fundamental axes of the Universe, similar to the matched circles 
searching in the circles-in-the-sky method. 

For practical application of our method, more realistic sim- 
ulations are necessary, and will be carried out in the next paper 



10 



H. Fujii & Y. Yoshii: An improved crystallography method for flat spaces 
Table 4. The results for different locations of the observer at {X, Y, Z). 



Location of 


Number of 


Number of 


Topological 


Topological 


observer 


quasars 


ghosts 


index 5 mult* 


index 5 simp* 


(0, 0, 0) 


2044 


24 


1.60 X 10-1° 


1.44 X 10-1" 


(0, 0.5, 0) 


2014 


18 


1.50 X 10-1° 


1.46 X 10-1" 


(0, 1,0) 


1970 


10 


1.46 X 10-1° 


1.47 X 10-1" 


(1,0, 0) 


1972 


24 


1.63 X 10-1° 


1.30 X 10-1" 


(3, 0, 0) 


2002 


26 


1.69 X 10-1° 


1.46 X 10-1" 


(5, 0, 0) 


1970 


26 


1.59 X 10-1° 


1.49 X 10-1" 



■5m„it and 5sin,p represent the topological index defined in equation (3) for the Hantzsche-Wendt space and for the simply connected Euclidean 

space, respectively. 



They will include more realistic characteristics of quasars such 
as spatial correlations, activity cycles, and anisotropic morphol- 
ogy. Technical problems such as magnitude limits and selection 
biases will also be considered there. One possible selection ef- 
fect relevant to our method, at least for translation, is clustering, 
which mimics a topological signal as was discussed in Marecki 
et al. (2005). They constructed the Lseiec filter, in the language 
of this paper, which removes a quadruplet [(Xi,Xj), {Xk,xi)] that 
is too compacted, i.e. satisfying \Xi - Xj\, \Xk - Xi\ < Lseiec- This 
filter will also be useful for our method. 

Presently available data in the latest versions of the Veron- 
Cetty & Veron quasar catalog (Velon-Cetty 8c Veron 2010) and 
the Sloan Digital Sky Survey (SDSS) quasar catalog (Schneider 
et al. 2010) will be used to make more precise constraints, when 
compared with the previous constraints that ignore the lifetime 
of quasars. Moreover, future observations that will detect hun- 
dreds of quasars with z> 6 (e.g., the Joint Astrophysics Nascent 
Universe Satellite (JANUS)), will enable us to remove the dis- 
agreement in the observational constraints using CMB data; 
specifically, we will detect or exclude the cubic 3-torus topol- 
ogy with L ^ 3.8L// mentioned by Aurich (2008). 

Acknowledgements. We thank Jeffrey Weeks and Adam Weeks Marano for 
beautiful figures of 17 multiconnected flat spaces. We also thank T. Minezaki, 
T. Tsujimoto, T. Yamagata, Y. Sakata, T. Kakehata, and K. Hattori for useful 
discussions and suggestions. 

Appendix A: Simulations for all types of 
holonomies 

We present the results for all types of holonomies here. The spa- 
tial resolutions used in these simulations are given in Table lA.ll 
b' is defined as the maximum value of s^, Sy, and e,, and we 
conservatively substitute (e', s', s') for s. Those for glide reflec- 
tion (type II) were determined in the same way as for half-turn 
corkscrew motion (type II), while the others were determined in 
the same way as for parallel translation. The reason we did so 
is that, for the latter types of holonomies, the deviation of the 
observer's location from the center always pushes all the faces 
away, and makes it harder to detect them. We cannot observe 
any ghosts except for those that stand very close to the faces of 
our Dirichlet domain. Hence a relative cosmic time between a 
quasar P,- and its ghost yiPi) is very small, which leads to small 
positional uncertainties (integrated effect) and small spatial res- 
olutions. 

We located the observer at the center and used the correct 
coordinate axes. The situations in which these quantities deviate 
from them are discussed in section 5.2. The results are given in 
Table lA.2l and Figures lATI and lATSl Relatively strong signals for 
half-turn corkscrew motion (type II) and glide reflection (type 



Table A.l. The spatial resolutions for each type of holonomies 
used here. 



Holonomy 


s (Mpc) 


e' (Mpc) 


Translation 


3.7 


6.0 


Half-turn corkscrew motion : I 


3.7 


6.0 


Half-turn corkscrew motion: II 


5.8 


7.3 


Quarter-turn corkscrew motion 


3.7 


6.0 


Third-turn corkscrew motion 


3.7 


6.0 


Sixth-turn corkscrew motion 


3.7 


6.0 


Glide reflection : I 


3.7 


6.0 


Glide reflection : II 


5.3 


6.5 



II) stem from the large spatial resolution. The signals for «-th 
corkscrew motion for n - 4, 3, and 6 are also strong, since we 
have to use two filters to detect y and y-i. 

We can see that the hills always appear, which are constituted 
by topological ghosts. For a case where the clear distinction be- 
tween the multiconnected space and the simply connected one 
is not seen, .s-,-histograms are indispensable to distinguish them. 
Our method is valid for all flat spaces. As mentioned in section 
5.1, zigzag features seen in the hills are not notable. They are 
merely stochastic and different from calculation to calculation. 
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Table A.2. The results for each type of holonomies. 



T-Tolonornv fvnp 


1 TTi hp r of 


1"! UillUCi KJl 




TnnnlncTipal 








iiiucA o niult 


itiHpy ^ ■ 


Translation : 1 


1980 


12 


6.60 X 10"" 


6.24 X 10"" 


Translation : 11 


1964 


28 


9.01 X 10"" 


6.07 X 10"" 


Translation ! Ill 


1970 


40 


9.59 X 10"" 


6.95 X 10"'° 


Half-turn corkscrew motion : 1 


1976 


16 


6.97 X 10"" 


7.21 X 10"" 


Half-turn corkscrew motion : 11 


2044 


24 


1.60 X 10"'o 


1.44 X 10"'° 


Quarter-turn corkscrew motion 


2007 


20 


1.36 X 10"'" 


1.43 X 10"'° 


Third-turn corkscrew motion 


2028 


12 


1.57 X 10"'" 


1.27 X 10"'° 


Sixth-turn corkscrew motion 


2014 


14 


1.42 X 10"'" 


1.37 X 10"'° 


Glide reflection : 1 


1973 


20 


7.76 X 10"" 


6.94 X 10"" 


Glide reflection : 11 


1950 


38 


1.34 X 10"'" 


1.11 X 10"'° 





Fig. A.l. The i, -histograms for each type of holonomy. Top left: 
one-pair translation; top right: two-pair translation; middle left: 
three-pair translation; middle right: half-turn corkscrew motion 
(type I); lower left: half-turn corkscrew motion (type II); lower 
right: quarter-turn corkscrew motion. 



Fig.A.2. Same as in Figure IaTI but for different holonomies. 
Top left: third-turn corkscrew motion; top right: sixth-turn 
corkscrew motion; lower left: glide reflection (type I); lower 
right: glide reflection (type II). 
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